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Abstract 

Nowadays various estimates of Lyapunov dimension of Lorenz-like systems 
attractors [1] are actively developed. Within the frame of this study the 
question arises whether it is possible to obtain the corresponding estimates 
of dimension for the Chen and Lu systems using the reduction of them to 
the generalized Lorenz system [2, 3]. In the work [4] Leonov’s method was 
applied for the estimation of Lyapunov dimension, and as a consequence 
the Lyapunov dimension of attractors of the Chen and Lu systems with the 
classical parameters was estimated. 

In the present work an inaccuracy in [4] is corrected and it is shown that 
the revised domain of parameters, where the estimate of Lyapunov dimension 
is valid, does not involve the classical parameters of the Chen and Lu systems. 
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1. Introduction 


Nowadays various estimates of Lyapunov dimension of attractors of gen¬ 
eralized Lorenz systems [1] 

{ X = a{y — x) 

y = cx + dy — xz (1) 

i = —bz + xy. 

are actively developed. In the work [4] the following result on the estimate 
of Lyapunov dimension of the generalized Lorenz system (1) is formulated. 

Theorem 1 ([4]). Let K he invariant compaet set of the generalized Lorenz 
system (1). Suppose that a > 0, b > 0, c + d > 0, {a — b){b + d) + ac > 0, 
and q = min(6, d). 

Then if it is satisfied one of the following conditions: 

1. d<0, a^c (4 - + 2a{b + d){2a - 3b) - b{b + d)^ > 0, 

2. d > 0, a > d, 

3. d = 0; and also one of the conditions: 

(a) 2a > b, af{3h + c) — b^ — ab{6b -1- c) > 0, 

(b) 2a < b, 3a?{h + c) — 6a6^ — b^ > 0] 

then dim^, K <3 - 2 {a+b-d) ^ valid. 

a—d+A/ {a+d)^-\-Aac 

However in the proof of this result in [4] there is an inaccuracy, the con¬ 
sideration of which the present work is devoted. 

2. Lyapunov dimension for the generalized Lorenz system 

Introduce a notion of Lyapunov dimension. Consider, for this purpose, 
the linearization of dynamical system (1) along a solution x(f, xq) 

^ = J(x(t,xo))u, (2) 

where J(x(t, xq)) is the (nxn) Jacobian matrix evaluated along the trajectory 
x(f,xo) of system (1). The fundamental matrix X{t,xo) of linearized system 
(2) is dehned by the variational equation 

X(t,xo) = J(x(t,xo))X(t,xo). (3) 
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Let (Ti(X(t,xo)) > ••• > (J„(X(t,xo)) > 0 be singular values of a fun¬ 
damental matrix X{t,xo) (the square roots of the eigenvalues of the matrix 
X{t,xo)*X{t,xo) are ordered for each t). 

Definition 1. The Lyapunov exponents (LEs/ in the point xq are the num¬ 
bers (or symbols ±cxDj 

LEi(xo) = limsup ^ lnai(X(f,xo)). (4) 

t^OO t 

Following [9], introduce a notion of Lyapunov dimension. 

Definition 2. A local Lyapunov dimension of a point xq in the phase space 
of dynamical system, generated by equation (1), is a number 


dim^xo = j(xo) 


LEi(xo) LEj(xo) 

|LEj+i(xo)| 


where LEi(xo) > > LEj(xo), are Lyapunov exponents; JIxq) G [l,n] 

the smallest natural number m such that 


( 5 ) 

is 


LEi(xo) LEm+i(xo) < 0, 


LEi(xo) LEm(xo) 

|LEm+i(xo)| 


< 1 . 


The Lyapunov dimension of invariant set K of dynamical system is de- 
hned by 

(LutilK = sup dim^XQ. (6) 


xqGK 


Note that LEs and Lyapunov dimension are invariant under the linear 
change of variable (see, e.g., [7]). 


^Two well-known definitions of Lyapunov exponents, which are used, are the following: 
the upper bounds of the exponential growth rate of the norms of linearized system solutions 
(LCEs) [5] and the upper bounds of the exponential growth rate of the singular values of 
fundamental matrix of linearized system (LEs) [6]. These definitions usually give the same 
values, but for a particular system, LCEs and LEs may be different (see, e.g., corresponding 
example in [7]). In general, it may be shown that the ordered LCEs majorise the ordered 
LEs, thus the values of Kaplan-Yorke (Lyapunov) dimension based on LEs and LCEs may 
be different. See, e.g., a more detailed discussion in [7, 8]. 

Note, that LEs are independent of the choice [7] of fundamental matrix at the point xq 
unlike the Lyapunov characteristic exponents (to get all possible values of LCEs, a normal 
fundamental matrix has to be considered [5]). 
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Alongside with widely used numerical methods for estimating and com¬ 
puting the Lyapunov dimension there exists analytical approach suggested 
by G.Leonov [10, 11, 12, 13, 1] that is based on the direct Lyapunov method. 

Let us outline the main theorems describing this approach . Let Ai (x, S) ^ 
■ ■ ■ ^ An(x, S) be the eigenvalues of the following matrix 

1 (SJ(x)S-‘ + (SJ(x)S-‘)*) . 

Theorem 2 ([14, 1]). Given an integer j G [1,?7.] and s G [0,1], suppose 
that there is a continuously differentiable scalar function d : M” —)■ M and a 
nonsingular matrix S such that 

Xi(x,S)-\ - \-Xj(x,S) + sXj+i(x,S)+ 'd(x) <0, VxGiL. (7) 

Then dim/, K ^ j + s. 

Here d is derivative of the function d with respect to the system (1). 

Theorem 3 ([15, 11, 13, 1]). Assume that there is a continuously differen¬ 
tiable scalar function d and a nonsingular matrix S such that 

Ai(x, S) + A 2 (x, S) + d(x) < 0, Vx G M’". (8) 

Then any solution of system (1), bounded on [0, -|-cxd), tends to an eguilibrium 
as t ^ -l-cxD. 


Thus, if condition (8) holds, then the global attractor of system (1) coin¬ 
cides with its stationary set. For proving Theorem 1 in [4] it is used above 
described approach with 


where 


and the function 
'd{x,y,z) 


S 


0 0 \ 
1 0 

\ 0 0 1 / 


p 


a 

{a — b){b d) ac 


(1 - s)V{x,y,z) ^ (1 - s)V{x,y,z) 

2^{a-2b-dr+i^ 2P{a + tiy-+ 4ac' 


(9) 


( 10 ) 


( 11 ) 
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where 


V{x,y,z) 



ax^zj +'y 2 y‘^ + {a^'yi + l 2 )z‘^ + 'y 3 xy + 'y 4 x‘^ -A^z. ( 12 ) 


The parameters 71 , 72 , 73,74 are chosen in dependence of system parameters 
in such a way that the conditions of Theorem 2 or 3 are satished. 

In proving Theorem 1 in the work [4] there is committed an inaccuracy, 
which does not permit one to obtain a condition for estimating the attractor 
dimension in the case d > 0. Namely, in the proof there is an inaccurate 
assertion that in the case 2 a + 6 — 2 d > 0 the inequality 


f a(2a + b — 2 d)^ 




(13) 


yields the inequality 

7 ^" < iV7l < 7f’, (14) 

where 7 , > 0, IV < 2a<i. 7 ^= {2a^ + ab)^i ^ 

Let us show that this assertion may be not true. Consider specific values 
of system parameters and running parameters for which the assertion is not 
satished. 

Take a = 1, 6 = 2, c = 6 , d = 1 for which 2a + b — 2d = 2, and 71 = 5, 
72 = —-jI. Then one has N <2. 

Condition (13) is satished since However 73 ^^ = 20 and 

consequently for N < 2 the relation 20 < 5N is impossible, i.e. condition 
(14) under such parameters is not satished. 

Correcting the above-mentioned inaccuracy, we obtain the following re¬ 
sult. 


Theorem 4. Denote by K the bounded invariant set of system (1), involving 
the point x = y = z = 0. Suppose that a > 0, 6 > 0, d > 0, c -|- d > 0, 
ac -|- (a — h)(h -|- d) > 0 , and there exists 7 such that the inequality 

—d{2a + b + 7 )^ -l- 8abd + A'yab > 0 


is satisfied. 
Let 


4b{—d{2a + b + 7 )^ -l- 8abd -|- 470 ^) 


p^{a + b + df 


46 


+ 


4a^(c-|-d) b{c + d) 


+ 
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+/9^(6 + d)(2a — 6 + 7 )^ 


< 0 , 


where p is defined as in (10). 
In this ease 

1. if 


a — d + 2b > 0, 
c<{b-d)(- + l], 


( 15 ) 


then any bounded on [0; +cxd) solution of system (1) tends to a eertain 
equilibrium as t ^ +cxd; 

2. zf 


then 


a — d -|- 5 > 0, 


0 {b - d) ^ + 1^ , 

(16) 

^ 2{a + b — d) 

a — d + + dfi + 4ac 

( 17 ) 


For proving Theorem 4 we uses the matrix S, dehned according to (9), 
and the function 

af N s)V{x,y,z) 

= 7:-TT-( 18 ) 


(a + dfi + 4ac]2 


where 


V{x,y,z) = 74a;^ + (-a7i+73)y^ + 73Z^ + — 7 ix'‘- 7 ia;^ 2 :- 7 i 72 a;?/-Y^. (19) 

4a b 

Exact value of dimension is obtained here from the comparison of dimen¬ 
sion estimate (8) and the value of Lyapunov dimension in zero point. The 
obtained value corresponds to the dimension of global attractor or a (local) 
B-attractor, involving zero equilibrium. 


3. Lyapunov dimension for the Chen and Lu systems 

In 2012 G.A. Leonov suggested to consider for the Chen and Lu systems 
the following change of variables [2] 

X —)■ hx, y —)■ hy, z hz, t h~^t (20) 
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Figure 1: Two-parametric domain for the Figure 2: Two-parametric domain for the 
Chen system Lu system 


with h = a. Under this transformation for a 7 ^ 0 the Chen and Ln systems 
are transformed to the generalized Lorenz system (1) with the parameters 

a —)■ 1 , 6 —d —( 21 ) 
a a 

Note that for a = 0 the Chen and the Ln systems become linear and their 
dynamics have minor interest. Thus, without loss of generality, one can 
assume that a = 1. The transformation (20) with h = a does not change 
the direction of time for the positive chaotic parameters considered in the 
works [16, 17]. Remark that the applying of transformations with the time 
inversion (see corresponding discussion of the Lorenz, Chen and Lu systems 
[3]) is not suitable for the study of Lyapunov dimension of invariant sets 
since the absolute values of Lyapunov exponents and corresponding local 
Lyapunov dimension in direct and backward time may be different (see, e.g., 
[ 8 , 18, 19, 3]). 

Below there are represented domains in the obtained two-dimensional 
parameter space (21) for which the condition of Theorem 4 is satisfied. On a 
plane there are also shown the points, corresponding to standard parameters 
of chaotic attractors. For the Lu system such parameters are (a, b, d) = 
(36,3,20), and after the change (b,d) = (A,h); for the Chen system — 
(a, b, d) = (35, 3, 28) and (&, d) = (^, |). 

As shown in Figs 1-2, nor the Lu system nor the Chen system standard 
parameters are situated in the obtained domain. While Lyapunov dimension 
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of the classical Lorenz attractor has been obtained analytically, for the Chen 
and Lu attractors it is still an open problem. Also by numerical simulation 
we do not hnd any chaotic attractor of the Chen and Lu systems with the 
parameters from the obtained domains. 

Applying the same approach, we extend the results, obtained in the case 
d = 0 in Theorem 1. 


Theorem 5. Suppose that a > 0,b > 0, c> 0, and the inequality ac + b{a — 
b) > 0 is satisfied. 

Let 7 ^^^^ > 0 be valid, where is the greater root of the equation 


p\2a-b + V^ + 16alr,{^^^ + L^±fi.-fifj-0, ( 22 ) 


where p is given according to (10) for d = 0. 
In this case 


1. if 


c <6 -+1 , 


(23) 


then any bounded on [0; +cx)) solution of system (1) tends to a certain 
equilibrium as t ^ +cxd; 

2. zf 

b 


c > 6 I-hi , 


then 


dimiK = 3 


2{a + b) 


+ + 4ac 


(24) 

(25) 


Consider, for example, the parameters a = l,b = 0.8, c = 9.2 for which 
there exists a chaotic attractor. For them 2a = 2, i.e. 2a > 6 is satished. 
However a^(36 + c) — b^ — ab{6b + c) = —0.112 < 0. Thus, for the above- 
mentioned parameters Theorem 1 cannot be applied. But for them Theorem 
5 can be applied since for them a higher root of equation (22) is equal to 
7^^^) ^ 8.583541432. 

Similarly it is possible to consider the generalized Lorenz system (1) not 
only with positive a, 6, but with the negative ones. 




Theorem 6. Suppose that K is a bounded invariant set of system 

{ X = a{y — x) 

y = rx — dy — xz (26) 

i = —bz + xy, 

involving the point x = y = z = 0. Suppose that a, 6 7 ^ 0, the inequality 
ra + {a — h){h — d) >0 is valid, and there exist numbers 71 , 72 , 73 such that 

71 > 0, (27) 



(28) 

(29) 

(30) 
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a + d + 2b > 0, 

ar < + ad + ab + bd, 


( 31 ) 


1 . if 


then any bounded on [ 0 ; + cxd ) solution of system ( 26 ) tends to a certain 
equilibrium as t ^ + cxd ; 

2 . zf 

a d b ^ 0, 

' a + d > 0, 
a + d < 0, 

a{r -d)>0, (32) 

' a + d + 2b < 0, 
a + d + 2b > 0, 
ar > b"^ + ad + ab + bd, 


then 


dimiK = 3 — 


2{a + b + d) 


a + d + A/(cr — dY + 4ar 


(33) 
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